Vertical cavity surface emitting terahertz laser 
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Vertical cavity surface emitting terahertz lasers can be realized in conventional semiconductor 
microcavities with embedded quantum wells in the strong coupling regime. The cavity is to be 
pumped optically at half the frequency of the 2p exciton state. Once a threshold population of 2p 
excitons is achieved, a stimulated terahertz transition populates the lower exciton-polariton branch, 
and the cavity starts emitting laser light both in the optical and terahertz ranges. The lasing 
threshold is sensitive to the statistics of photons of the pumping light. 



PACS numbers: 78.67.Pt,78.66.Fd,78.45.+h 

Creation of efficient sources of terahertz (THz) radia- 
tion is of high importance for various fields of modern 
technology, including information transfer, biosensing, 
security and others 1 . These applications are currently 
limited due to the lack of compact and reliable solid state 
sources of THz radiation. The main obstacle preventing 
the creation of such a source is the low rate of spon- 
taneous emission of THz photons: according to Fermi's 
golden rule this rate is proportional to the cube of the fre- 
quency and for THz transitions should be roughly tens of 
inverse milliseconds, while lifetimes of crystal excitations 
typically lie in the picosecond ranged. The attempted 
strategies to improve this ratio include using the Purcell 
effect 4,5 by embedding the sample inside a THz cavity, 
or using the cascade effect in quantum cascade lasers^ 
(QCLs). Nevertheless, until now QCLs in the spectral 
region around 1 THz remain costly, have macroscopic 
dimensions, are short-lived and still show a quantum ef- 
ficiency of less than 1%. 

Recently, it has been proposed that the emission rate 
of THz photons may be increased by bosonic stim- 
ulation if the THz transition feeds a condensate of 
exciton-polaritons. In Refs0-[il the authors suggest us- 
ing the transition between the upper and lower polariton 
branches in a semiconductor microcavity in the polariton 
lasing regime. The radiative trasition between two polari- 
ton modes, accompanied by emission of a THz photon, 
may become allowed if an electric field mixing polariton 
and dark exciton states is applied to the cavity. THz 
photons would be emitted in the plane of the cavity, and 
a supplementary lateral THz cavity would be needed to 
provide the positive feedback. 

Here we propose an alternative model of a microcavity 
based THz laser, which has significant advantages com- 
pared to one based on the transition between two po- 
lariton branches. We propose using two-photon pump- 



ing of a 2p exciton state, as has been realised already 
in GaAs based quantum well structure s 10 ^ 1 . The direct 
transition to or from the 2p exciton state with emission 
or absorption of a single photon is forbidden by optical 
selection rules. Instead, a 2p exciton can radiatively de- 
cay to the lower exciton-polariton mode formed by the 
Is exciton and cavity photon. This transition is accom- 
panied by emission of a THz photon^. The inverse pro- 
cess (THz absorption by a lower polariton mode with 
excitation of a 2p exciton) has been recently observed 
experimentally^. The THz transition from the 2p state 
pumps the lowest energy exciton-polariton state, which 
eventually leads to the polariton lasing effect, widely dis- 
cussed in the literature 1 ^. A macroscopic occupation of 
the lowest energy polariton state stimulates emission of 
THz photons, so that, in the polariton lasing regime, the 
cavity would emit one THz photon for each optical pho- 
ton emitted by the polariton laser, ideally. 

The design of the laser and the involved energy levels 
are illustrated schematically in Figure 1. This design has 
two crucial advantages with respect to that previously 
considered 7 : it allows for operation with an optically al- 
lowed THz transition, and it provides vertical emission 
of THz photons. The whole structure is microscopic; 
no waveguides or THz cavities are needed. This latter 
point is also a significant advantage with respect to the 
quantum cascade laser, which operates in the wave-guide 
geometry. Moreover, the suggested scheme is very inter- 
esting from a fundamental point of view, since, as shown 
below, the threshold of the proposed laser appears to be 
sensitive to the statistics of photons of the pumping light. 
Here we present the quantum model of this laser based 
on the Liouville equation for the density matrix describ- 
ing the 2p exciton state, the lower polariton mode, and 
optical and THz photons. 

The system under study consists of 2p excitons and Is 
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In the particular system which we consider, the coherent 

and incoherent parts of the Hamiltonian can be written 
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FIG. 1: (a) Schematic of the polariton dispersion relation, 
showing the lower-polariton (LP) and upper-polariton (UP) 
branches, as well as the 2p exciton state with frequency uj p . 
The pump frequency, cj a , is half that of the 2p exciton state, 
as discussed in the text, (b) The structure considered in this 
work: a semiconductor microcavity consisting of an active 
layer containing quantum wells (QWs) sandwiched between 
two distributed Bragg reflectors (DBRs). The structure is 
pumped vertically, i.e. in the direction perpendicular to the 
microcavity plane, and the resulting THz emission from the 
cavity is in the same direction. 



exciton-polaritons interacting with an external electro- 
magnetic field and phonons in a semiconductor microcav- 
ity. Hence, our consideration involves both coherent and 
incoherent processes. Consequently, it is convenient to 
split the total Hamiltonian of the system into two parts, 



H — H r + H d 



(i) 



where the part H c includes all types of coherent pro- 
cesses in our system and part Hd describes the deco- 
herence in the system due to exciton interactions with 
acoustic phonons, modelled as a classical reservoi r 15 ! 16 . 

The equation for the density matrix can be written in 
the following form: 



dg i 



■Lq, 



where L is a Lindblad superoperator defined below. The 
incoherent part of the Hamiltonian can be divided into 
two parts, Hd = + H~), where H~^~ creates an 

excitation in the quantum system (and thus annihilates 
an excitation in the classical reservoir), and H~ , con- 
versely, annihilates an excitation in the quantum system 
and creates an excitation in the classical reservoir. The 
subscript index i numerates different types of decoher- 
ence processes in the system. 

The Lindblad terms can be then written a s 16, 17 
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J2{(H+pH-+H-gH+) 



(2) 



(H+H- + H-H+) g - g (H+H- + H-H+) } , 



where S(AE) accounts for energy conservation, and in re- 
alistic calculations should be taken as an average inverse 
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Here p and s denote annihilation operators of the 2p ex- 
citon and the lowest energy Is polariton states, respec- 
tively, a is the annihilation operator for laser photons ex- 
citing the 2p exciton state, and c is an annihilation opera- 
tor for THz photons produced by the 2p —> Is transition. 
e p , e s , e a and er are the energies of the 2p exciton, Is 
exciton-polariton, pump photon (2e a = e p ) and THz pho- 
ton, €t = e p — e s , respectively. Constants g and G define 
the strengths of the two-photon absorption and Is — >• 2p 
radiative transitions, respectively. While the Is — >• 2p 
transition strength is determined by the dipole matrix 
elements, the main contribution to the two-photon 2p 
excitation comes from processes with intermediate states 
in the valence or conduction bands 18 . We assume that 
the upper polariton state is higher in energy than the 2p 
exciton state, and thus can be excluded from consider- 
ation. The classical reservoir consists of the photons of 
the external laser light and THz photons. Note, that the 
coherent part of our Hamiltonian describes only the free 
modes of the system and thus is irrelevant to the kinetic 
equations. 

Now, for the occupancy of the 2p exciton mode one 
can obtain after straightforward algebra: 
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~dt~ 
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Re{Tr(^-;[p + p;i7 + ]])} 



= W a 



i(a +2 a 2 )-(p+p(2a+a + l)) 



+W G [(s + sc + c) - (p+p(s+s + c + c + 1))] 



(9) 



where (...) = Tr {...£>} denotes averaging with the ap- 
propriate density matrix g, W g = 4g 2 /h( and Wg = 
2G 2 /h(. The equations for the polariton mode occu- 
pancy N s = and the terahertz mode occupancy 
N c = (c + c) are analogous to those for N p . The occu- 
pancy of the pumping mode N a = (a + a) is defined by 
the intensity of the external pump and we do not need 
to write an independent dynamic equation for it. The 
same holds true for the higher order correlators involv- 
ing pump operators, e.g. (a +2 a 2 ). It follows from Eq. (J9j) 
that the dynamic equations for the occupancies of the 
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modes contain quantum correlators of fourth order, such 
as (s+sc+c). For them one can also write the dynamic 
equations analogous to Eq. ((9]), which would contain cor- 
relators of sixth order. Proceeding further, one would 
obtain an infinite chain of coupled equations for the hi- 
erarchy of correlators. 

In order to solve this chain of equations one needs to 
truncate the correlators at some stage. Here we use the 
mean-field approximation, which consists in truncation of 
the fourth-order correlators into products of second-order 
ones. One can approximate (p+pa+a) « (p+p)(a+a) = 
N p N a , (s+sc+c) w N S N C , ... etc. 

A particular analysis is needed for the truncation of 
the correlator (a +2 a 2 ) containing four operators corre- 
sponding to the pumping mode. Using the definition of 
the second order coherence g^ 2 \0), it can be represented 
as 
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(a+V)=^ 2 )(0)iV a 2 . 



(10) 



Finally, the closed set of equations of motion describing 
the dynamics of our system reads 



dN p 
~dt~ 

dN s 
dt 



dNc 
dt 



N 2 a - N p (2N a + 1) 



(11) 



+W G {N S N C (N p + 1)-N P (N s + 1) (N c + 1)} , 
N 

-— ~ (12) 

-W G {N s N c (N p + 1)-N p (N s + 1) (N c + 1)} , 
A 

-— - (13) 
-W G {N s N c (N p + 1)-N p (N s + 1) (N c + 1)} , 



where we have introduced the non-radiative lifetime of 
the 2p exciton state , r p , and lifetimes of lower polaritons 
and THz photons, r s and r c , respectively. 

Let us analyze in more detail the kinetic equation for 
the 2p state occupation N p . The terms describing incom- 
ing and outcoming scattering rates for this state can be 
written as 

(l + N p )(a +2 a 2 )-N p (a 2 a +2 ). 

As we already mentioned in Eq. ([TU]h the incoming rate 
is proportional to g^ 2 \0)N 2 . The outgoing rate can be 
recast as oc (a 2 (a+) 2 ) = g( 2 \0)N 2 + AN a + 2. Hence, 
all terms containing N 2 N P cancel each other and Eq. 
(fTTj) holds. Note, that for processes of non-degenerate 
two-photon absorption or emission neglected here, two- 
photon emission processes where photon frequencies are 
different can be included in r p , the decay rate of the 
2p state. The scattering rates have standard form (1 + 
N p )N ai N a2 - N p (N ai + l)(iV a2 + 1), independent of the 
photon field statistics^. In what follows we assume that 
N a ^> 1 making it possible to neglect all other modes. 

The key feature of Eq. ([TT]) is that the pumping term 
W g g( 2 \0)N% contains the second order coherence of the 
pump, i.e. strongly depends on its statistics. For a laser 
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FIG. 2: Occupation of the p- and s-states in the steady state 
(left axis) and terahertz generation rate (right axis), as a func- 
tion of pumping mode occupancy. The subscript 1 or 2 refers 
to the value of # (2) (0) taken. W g = 1, W G = 10, l/r s = 



1000, 1/t p = 5000, in units of W g 



pump corresponding to the coherent state this statistics 
is Poissonian and ^^ 2 ^(0) = 1. On the other hand dif- 
ferent radiation sources may provide different statistics 
of pump photons and different values of g( 2 )(0). For ex- 
ample, for a thermal pump g( 2 )(0) = 2. Different val- 
ues of g( 2 )(0) lead to the different lasing thresholds in 
the system. Note that dependence of two-photon opti- 
cal processes on the second order coherence of light has 
been discussed by Ivchenko in Ref. |2(J In Figure 2, we 
plot the steady state solutions for N p and N s , assuming 
that the terahertz mode occupation is zero, for the 
cases of coherent and thermal pumps. For a qualitative 
understanding, we note that the threshold is reached if 
N s ~ 1 5 or equivalently when N p ~ 1/t s Wg, since this 
signals the onset of Bose stimulation of the transition 
from the 2p-state to the lower polariton, accompanied by 
terahertz emission. The terahertz generation rate is given 
by T = W G N P (N S + 1) « l/r s A^ s , and therefore shows 
the same threshold behaviour with pumping intensity as 
Ag, as can also be seen in Figure 2. We observe that the 
threshold is higher in the case of a coherent pump, which 
may be understood in terms of the enhanced losses from 
the 2p-state due to stimulated two-photon emission. 

The dynamics of switching the system to polariton and 
terahertz lasing is essentially governed by the dynamics of 
the correlator (a +2 a 2 ), i.e. by the product of squared in- 
tensity and second order coherence of the pumping light. 
The efficiency of two-photon absorption W g and the Is 
polariton lifetime r s , dependent on the quality factor of 
the optical cavity, are two main factors which govern the 
threshold to lasing. The polarisation of pumping defines 
the polarisation of emission, and the frequency of THz 
radiation may be tuned by changing the exciton-photon 
detuning in the microcavity. Interestingly, THz lasing in 
our system is possible in the absence of a THz cavity. In 



order to demonstrate this, we have taken the terahertz 
photon lifetime r c = in the calculation shown in Figure 
2. One can see that THz lasing begins simultaneously 
with polariton lasing. The stimulated emission of THz 
photons is achieved due to the bosonic stimulation of this 
transition by occupation of the final (Is polariton) state. 
The system emits THz radiation through the Bragg mir- 
rors of the microcavity, which are transparent at THz 
frequencies. 

In conclusion, we have developed a quantum theory 
of vertical cavity surface emitting terahertz lasers. The 
vertical terahertz lasing coexists with polariton lasing 
and, consequently, is characterised by a low pumping 
threshold. Interestingly, the value of the threshold is 
strongly dependent on the statistics of photons of the 
pumping light. For practical realization of compact ter- 
ahertz lasers operating at room temperature, microcav- 
ities based on wide band gap semiconductors (GaN or 
ZnO) would seem to be the most advantageous. The res- 



onant two photon pumping of 2p exciton states in such 
a cavity may be assured by a conventional vertical cav- 
ity light emitting diode (VCLED) emitting in red. This 
prospective structure would consist of a GaN or ZnO mi- 
crocavity grown on the top of a GaAs based VCLED 
structure. As emission of coherent terahertz light co- 
exists with the polariton lasing in the optical frequency 
range in this scheme, one can use the visible blue or green 
light produced by the polariton laser as a marker for the 
terahertz light beam, which may be important for appli- 
cations in medicine, security control, etc. 
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